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ABSTRACT 


The  reduction  of  test  time  in  low  pressure  shock  tubes!  due 
to  a  laminar  wall  boundary  layer,  has  been  analytically  in¬ 
vestigated.  In  previous  studies  by  Roshko  and  Hooker  the 
flow  was  considered  in  a  contact  surface  fixed  co-ordinate 
system.  In  the  present  study  it  was  assumed  that  the  shock 
moves  with  uniform  velocity,  and  the  flow  was  investigated 
in  a  shock  fixed  co-ordinate  system.  Unlike  the  previous 
studies,  the  variation  of  free  stream  conditions  between  the 
shock  and  contact  surface  was  taken  into  account.  It  was 
found  that  p,  a  parameter  defined  by  Roshko,  is  considerably 
larger  than  the  estimates  made  by  Roshko  and  Hooker  except 

for  very  strong  shocks.  Since  test  time  is  proportional  to 
-2 

p  ,  previous  estimates  of  test  time  are  too  large,  particu¬ 
larly  for  weak  shocks.  The  present  estimates  for  p  appear 
to  agree  with  existing  experimental  data  to  within  about  10 
percent  for  shock  Mach  numbers  greater  than  5.  In  other 
respects,  the  basic  theory  is  in  general  agreement  with  the 
previous  results  of  Roshko. 
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I.  INTRODUCTION 


In  an  ideal  shock  tube  (i.  e.  ,  neglecting  wall  and  real  gas  effects),  the  shock 
and  the  contact  surface  both  move  with  a  constant  velocity  and  the  flow  between 
them  is  uniform  (Fig.  1).  In  an  actual  shock  tube  flow,  however,  the  presence 
of  a  wall  boundary  layer  causes  the  shock  to  decelerate,  the  contact  surface 
to  accelerate,  and  the  flow  to  be  nonuniform  (e.  g.  ,  Ref.  1).  Analytical  and 
experimental  studies  of  the  nonuniformities  in  shock  tube  flows  have  been 
presented  in  References  1  through  10  (as  well  as  by  others). 

The  analytical  theory  presented  in  References  5  and  6  is  applicable  when  the 
wall  boundary  layer  introduces  only  small  perturbations  into  the  ideal  flow. 

In  these  references,  it  was  shown  that  the  wall  boundary  layer  between  the 
shock  and  the  contact  surface  acts  as  an  aerodynamic  sink,  removing  mass 
from  the  region  between  the  shock  and  the  contact  surface.  This  mass  removal 
causes  the  shock  to  decelerate  and  the  contact  surface  to  accelerate,  in  agree¬ 
ment  with  the  experimental  observations  of  Reference  1. 

As  the  length -to -diameter  ratio  of  a  shock  tube  is  increased  and  as  the  initial 
pressure  in  the  low  pressure  section  is  reduced,  the  wall  boundary  layer 
effects  become  more  pronounced.  (The  study  of  dissociation  and  ionization  in 
shock  tubes  has  stimulated  the  use  of  low  initial  pressures.  )  Duff,  ^  in  a  study 
of  the  flow  about  12  feet  from  the  diaphragm  in  a  1-1/8  in.  diameter  shock  tube 
with  initial  pressures  of  the  order  of  1  mm  Hg,  found  the  flow  to  be  strikingly 
different  from  that  in  a  conventional  shock  tube.  In  a  conventional  tube,  the 
separation  between  the  shock  and  the  contact  surface  (and  therefore  the  test 
time)  increases  with  distance  from  the  diaphragm.  However,  in  a  low  pressure 
shock  tube,  Duff  observed  that  the  separation  reaches  a  limiting  value  and 
remains  constant  with  distance  thereafter.  When  this  limiting  condition  is 
reached,  the  shock  and  the  contact  surface  both  move  with  equal  and  constant 
velocity.  This  phenomenon  must  be  taken  into  account  when  estimating 
test  time  in  low  density  shock  tubes. 
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(b)  FLOW  AT  TIME  i, 


Shock  Tube  Flow  in  Laboratory  Coordinates 


Duff  correctly  explained  the  limiting  flow,  where  shock  and  contact  surface 

move  with  equal  velocity,  as  one  wherein  the  flow  passing  through  the  shock 

equals  the  flow  which  moves  past  the  contact  surface  due  to  the  wall  boundary 

layer.  He  referred  to  the  contact  surface  as  a  "leaky  piston.  "  This  effect 

was  later  studied  experimentally  and  analytically  by  Roshko.  In  the  analytical 

portion  of  Reference  8,  Roshko  considered  the  flow  in  a  co-ordinate  system  in 

which  the  contact  surface  was  stationary.  The  separation  distance  between 

the  shock  and  contact  surface  at  any  instant  was  found  from  a  mass  balance 

that  equated  the  mass  flow  through  the  shock  to  the  sum  of  the  boundary  layer 

mass  flow  moving  past  the  contact  surface  and  the  rate  of  accumulation  of  mass 

9 

between  the  shock  and  contact  surface.  (Anderson  had  used  a  similar  approach.  ) 
Roshko  used  the  boundary  layer  theory  presented  in  References  11  and  1Z 
(which  applies  for  the  boundary  layer  behind  a  shock  moving  with  uniform 
velocity).  He  also  developed  an  approximate  boundary  layer  theory  to  account 
for  real  gas  effects.  (The  latter  are  treated  more  accurately  in  Refs.  13 
and  16.  )  Experimental  measurements  of  test  time  for  a  variety  of  initial 
pressures  were  obtained  by  Roshko,  which  confirmed  the  basic  theory. 


Hooker^  noted  an  erroneous  velocity  transformation  in  Reference  8  and  also 
noted  that  Roshko  had  left  out  a  term  in  the  mass  balance  equation,  namely,  the 
accumulation  of  mass  between  the  shock  and  contact  surface  due  to  nonuniform 
density  in  the  boundary  layer.  He  proceeded  to  correct  Roshko's  theory  for 
these  effects  and  claimed  a  somewhat  improved  correlation  between  theory  and 
experiment.  However,  the  improvement  was  not  dramatic. 


The  analytical  results  of  Roshko  and  Hooker  describe  the  basic  features  of  the 
flow  in  a  low  density  shock  tube  and  are  widely  used  to  estimate  test  time. 
However,  the  flow  model  used  in  both  papers  (i.  e.  ,  a  co-ordinate  system  in 
which  the  contact  surface  is  fixed)  contains  several  inherent  contradictions.  For 
example,  both  authors  solve  for  the  nonuniform  shock  velocity  relative  to  the 
contact  surface  but  assume,  for  one  term  in  the  mass  balance  equation,  that  the 
shock  velocity  is  constant.  The  nonuniformity  of  the  flow  between  the  shock  and 
the  contact  surface  is  not  taken  into  account  since  both  authors  use  a  boundary 
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layer  theory  based  on  uniform  flow  between  the  shock  and  the  contact  surface. 
They  recognized  that  the  latter  assumption  becomes  correct  only  in  the  case  of 
very  strong  shocks  but  made  no  attempt  to  modify  the  theory  for  the 
"not -so-strong"  shock  despite  the  fact  that  much  of  their  experimental  data  were 
obtained  at  moderate  Mach  numbers. 

In  the  present  paper,  the  problem  of  test  time  in  a  low  density  shock  tube  is 
investigated  by  considering  a  flow  model  wherein  the  shock  moves  with  uniform 
velocity.  The  co-ordinate  system  is  one  in  which  the  shock  is  stationary,  and 
the  wall  moves.  The  present  model  is  self-consistent.  A  boundary  layer 
theory  is  developed  to  take  into  account  the  nonuniform  flow  between  the  shock 
and  the  contact  surface.  The  solution  applies  to  shocks  of  moderate  strength, 
as  well  as  to  strong  shocks.  Previous  experimental  data  are  re-examined  in 
the  light  of  the  present  theory. 
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II.  STEADY  STATE  SOLUTION 

The  experimental  results  presented  in  References  7,  8,  and  10  indicate  that 
the  shocked  gas  in  a  low  density  shock  tube  ultimately  reaches  a  steady  state 
condition  where  both  the  shock  and  the  contact  surface  move  with  equal  and 
constant  velocity.  The  flow  between  the  shock  and  the  contact  surface  is 
then  steady  when  viewed  in  a  co-ordinate  system  in  which  the  shock  and  the 
contact  surface  are  stationary.  In  this  co-ordinate  system  the  wall  moves 
with  velocity  uw  (which  equals  the  shock  velocity  Ug  in  the  laboratory 
system).  This  steady  flow  is  investigated  herein  with  the  primary  object 
of  determining  the  separation  distance  between  the  shock  and  the  contact 
surface.  The  problem  of  unsteady  flow,  where  the  shock  and  the  contact 
surface  have  different  velocities,  is  treated  in  Section  III. 

Steady  flow  is  illustrated  in  Fig.  Z.  The  shock  is  located  at  I  =  0,  and  the 
free  stream  portion  of  the  contact  surface  at  I  =  1  The  flow  upstream  of 
the  shock  is  denoted  by  subscript  oo  and  moves  with  velocity  uw»  as  does 
the  wall.  Free  stream  conditions  between  the  shock  and  the  contact  surface 
are  denoted  by  subscript  e.  Free  stream  conditions  directly  downstream  of 
the  shock  have  the  additional  subscript  o.  The  percentage  of  mas 8  flow  in 
the  boundary  layer  increases  with  l  such  that  all  the  mass  flow  is  in  the 
boundary  layer  at  I  and  the  free  stream  is  stationary  at  that  location. 

g 

Roshko  obtained  an  estimate  for  f  in  the  following  manner.  The  flow 
rate  through  the  shock,  mg,  equals 

rh  -  (p  u  )  A  (la) 

s  e  e  o 

where  A  is  the  cross-sectional  area  of  the  tube.  Roshko  assumed  that  the 
boundary  layer  was  thin  and  characterized  the  flow  in  the  boundary  layer  at 
the  contact  surface,  mc,  by 
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Fig.  2.  Flow  Between  Shock  and  Contact  Surface 
In  Shock-Stationary  Coordinate  System 
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m  =  Lp  (u  -  u  )6d  (lb) 

c  rw,  o'  w  e,  o'  R  '  1 

where  L  is  the  perimeter  of  the  tube;  6^  is  a  characteristic  boundary  layer 

thickness  at  f  ;  and  p  and  u  -  u  are  characteristic  boundary  layer 
m  rw,  o  w  e,  o  7  7 

densities  and  velocities,  respectively.  The  boundary  layer  thickness  was 
further  characterized  by 


6 


R 


(1c) 


where  p  is  a  constant  that  must  be  found  from  an  exact  solution  of  the  bound¬ 
ary  layer  development  in  the  tube.  Equating  Eqs.  (la)  and  (lb),  since  the 
flow  is  steady,  yields  the  following  expression  for  I  ; 


f 

m 


u  u 

e, o  e, o 


u  -  U  V 

w  e,  o  w,  o 


(2) 


where  d  =  4A/L.  is  the  hydraulic  diameter.  Assume  that  the  temperature 

upstream  of  the  shock  is  at  a  standard  condition  so  that  T  =  T  .,  a  =  a  .  , 

and  p  =  Also,  assume  that  the  wall  remains  at  its  initial  temperature 

so  that  T  =  T  .  From  continuity,  p  u  =  (p  u  )  .  Equation  (2)  can  then 
w  s  t  oo  w  e  e  o 

be  put  in  the  form 


^st  *  m 
d2 


oo 


WM 


l6p 


2  P 


e,  o 


W 


-ll^st 


(3) 


where  W  =  u  /u  =  p  /p  ,  M  =  u  /a  ,  and  p  .  is  a  standard  pressure 
w  e,o  e ,  o  r oo  s  w  oo  st  r 

(usually  an  atmosphere).  The  right  hand  side  of  Eq.  (3)  depends  primarily 
on  Mg.  Hence,  for  a  given  Mgl  the  separation  distance  I  is  proportional 


-7- 


to  the  product  d  p  .  This  can  result  in  very  short  test  times  for  shock 
r  oo  7 

tubes  with  low  initial  pressure. 

Equation  (3)  does  not  yield  numerical  results  for  i  unless  an  accurate 

estimate  of  (3  is  available.  Estimates  of  p  have  been  presented  in  References  8, 

10  nnH  1  4  .  T  h  f*  nrimarv  rmrnr»<?o  of  tV\.»  nroqnnt  ic  to  i  mni*Avp  tKp  a 

•  i  —  ■  -  -  i-  ;  r  —  x~  —  -  —  * ' '  r  ~  —  ^  *  “ —  • —  *  ~  *jr  -  w  '  w  **  — 

estimates,  particularly  tor  flows  with  moderate  M  .  This  is  done  by  taking 

s 

Eq.  (2)  as  the  defining  relation  for  p 


f  1 6f 


e  i  o 


u 


m  V  w,  o> 


W  -  1 


e,  o 


w,  o 


(4) 


and  finding  i  as  accurately  as  possible  from  a  consideration  of  the  bound¬ 
ary  layer  development  in  the  flow  illustrated  in  Fig.  2.  A  first  estimate  is 
made  below  by  considering  the  boundary  layer  to  develop  in  a  uniform  external 
stream.  An  improved  estimate  is  then  made  by  employing  the  concept  of 
local  similarity. 

A.  UNIFORM  EXTERNAL  FREE  STREAM  APPROXIMATION 

Boundary  layer  development  for  the  case  of  an  external  free  stream  that  does 

not  vary  with  l  is  illustrated  in  Fig.  3  and  is  discussed  in  the  Appendix.  Let 

l  correspond  to  the  value  of  i  at  which  the  excess  mass  flow  in  the  boundary 

layer  equals  the  mass  flow  entering  through  the  shock.  This  gives 

A(p  u  )  =  L(p  u  )  (-6*)  where  6*  is  the  boundary  layer  displacement 

e  e  o  e  e  (3 

thickness  at  I  .  Combining  the  latter  expression  with  Eqs.  (4)  and  (A-4) 
(Appendix)  yields 


The  subscript  o  has  been  added  to  p  to  indicate  that  this  value  is  based 
on  a  bounday  layer  with  a  uniform  external  stream. 

> 

f 
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Fig.  3.  Boundary  Layer  Development  Assuming  Uniform  Free  Stream 


-9- 


Numerical  values  of  Pq  have  been  computed  assuming  an  ideal  gas  and 
pp  =  constant  in  the  boundary  layer  (Appendix).  These  are  listed  in  Table  1 
for  Prandtl  number  (<r)  =  0.  72  and  1;  y  =  7/5  and  5/3;  and  for  various  shock 
strengths. 

A  solution  for  the  boundary  layer  behind  a  strong  shock  moving  with  uniform 
velocity  into  air  at  an  initial  pressure  of  0.  001  atmospheres  and 
T  =  T  =  522°R  has  been  presented  in  Reference  13.  (Some  results  at  0.  01 

W  OO 

atmospheres  were  also  obtained.)  This  is  a  "real  gas"  solution  in  the  sense 
that  equilibrium  gas  charts  were  used  to  obtain  the  flow  across  the  shock  and 
in  the  boundary  layer.  The  Sutherland  relationship  was  used  to  evaluate  the 
variation  of  pp.  in  the  boundary  layer.  Some  of  the  boundary  layer  results 
are  given  in  Table  2.  The  corresponding  values  of  pQ  are  included  in  this 
Table. 

It  is  expected  that  these  values  of  pQ  will  overestimate  particularly  for 

shock  Mach  numbers  that  are  not  large.  This  is  due  to  the  fact  that  the 

relative  velocity  between  the  wall  and  the  free  stream  increases  from 

u  -  u  at  I  =  0  to  u  at  I  (compare  Figs.  2  and  3).  Hence  the  excess 

mass  flow  in  the  boundary  layer  will  be  greater  at  a  given  l  than  that  obtained 

from  the  above  model  (which  assumes  the  relative  velocity  to  remain  constant 

at  u  -  u  ).  This  will  result  in  smaller  i  and  larger  p  than  obtained 
w  e,o'  m  ©  r 

from  Eq.  (5a).  However,  for  very  strong  shocks,  where  u  is  small 

e ,  o 

relative  to  uw>  Eq.  (5a)  should  give  accurate  results. 

In  Reference  8,  Roshko  used  two  different  expressions  to  evaluate  p.  The 
first  expression  arose  from  a  velocity  transformation  error  (pointed  out  by 
Hooker)  and  is  the  same  as  Eq.  (5a).  His  results  for  cr  =  1,  y  =  7/5  (in  Table  I 
of  Ref.  8)  are  in  agreement  with  Table  1  herein.  He  also  used  the  expression 

E  >/w  -  1  o 
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which  is  the  basis  for  the  results  he  presents  in  his  Tables  II  to  IV.  Since 
I^/ff  -  ri)^  is  positive,  p^  is  somewhat  smaller  than  pQ  and  therefore 
leads  to  greater  error  when  computing  I  ,  Equation  (5b)  is  the  expression 
that  arises  naturally  from  the  contact  surface  fixed  co-ordinate  system 
used  in  Reference  8  and  is  presumably  the  expression  Roshko  meant  to  employ 
throughout  his  paper.  In  References  10  and  14,  Eq.  (5b)  was  therefore  used 
to  evaluate  p.  This  is  discussed  further  in  Section  IV. 

B.  LOCAL  SIMILARITY  APPROXIMATION 

In  the  present  section,  the  streamwise  variation  of  free  stream  properties 
due  to  the  increase  in  boundary  layer  mass  flow  with  I  is  taken  into  account. 
The  development  of  the  boundary  layer  and  the  variation  in  free  stream 
properties  are  treated  simultaneously.  The  boundary  layer  growth  is  found 
by  assuming  that  at  each  station  it  is  similar  to  a  corresponding  boundary 
layer  developing  in  a  uniform  free  stream  behind  a  shock  moving  with  uniform 
velocity  (i.  e.  ,  local  similarity). 

Since  the  flow  is  steady  (Fig.  2),  the  net  mass  flow  through  the  shock  equals 
the  net  mass  flow  at  any  station  f .  Thus 


f  <» 

A(p  u)  =  Ap  u  +  L  I  (pu-p  u  )  dy 
xre  e'o  re  e  Ja  e  e'  7 


(6) 


In  Eq.  (6),  it  is  assumed  that  the  boundary  layer  thickness  is  small  compared 
with  d;  thus  the  integrand  is  nonzero  only  in  the  region  close  to  the  wall. 
Define 


6 


.  p  u 
4  e  e 

d (peue)Q 


(-6*) 


(7a) 
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where  6*  is  the  local  boundary  layer  thickness 


s*  ■/. 

Note  that  6  is  the  ratio  of  the  excess  mass  flow  through  the  boundary  layer 
at  i,  to  the  mass  flow  through  the  shock.  Thus  6  varies  from  0  at  i  =  0 
to  1  at  i  =  Equations  (6)  and  (7)  then  give 


6  = 


P  u 
re  e 

(peue)0 


(8) 


which  relates  the  free  stream  conditions  to  the  local  boundary  layer  displace¬ 
ment  thickness. 

The  concept  of  local  similarity  is  now  introduced.  It  is  assumed  that  the 
boundary  layer  profile  at  each  f  corresponds  to  that  for  a  boundary  layer 
associated  with  a  uniform  free  stream  (equal  to  the  local  free  stream)  and 
a  wall  velocity  u^,  The  origin  of  this  boundary  layer  is  at  I which  is 
initially  an  unknown  function  of  I.  (See  Fig.  4.)  The  origin  I  ^  is  chosen 
such  that  the  excess  flow  in  the  boundary  layer  at  each  I  has  the  correct 
local  value.  It  is  also  assumed  that  the  boundary  layer  growth  at  each 
section  is  the  same  as  that  for  the  corresponding  uniform  free  stream 
boundary  layer. 

The  local  displacement  thickness  is  then  (Eq.  A-4) 


/2t< 


i 


+  1  ] 

oo 


(9) 
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4 

4.  Boundary  Layer  Development  With  Uniform  Free 

Stream,  Wall  Velocity  u  ,  and  Origin  at  I. 

w  i 
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where  (f  -  n)  and  I  are  functions  of  the  local  free  stream  and  wall 

’  OO  00  _ 

conditions  (Appendix  ).  From  Eq.  (7a),  6  becomes 


where 


,syR“. 


(10a) 


e. life)  I  isi 


fdVpe/oJ  Ue,  o 


(10b) 


:  2 


f  4  /Pw\  1  vw,  o  | 

Ld\Pe/oJ  ue,  o  1 


(10c) 


/  Pe 

H  ^./V— — 

e  V  pe.0 


[(f-Tl)  +1  ] 

L'  ''oo  00  J 


(lOd) 


V  s  u  /u 

e  e ,  o 


(10e) 


The  nondimens ional  variable  £  now  replaces  i.  In  deriving  Eq.  (10a)  from 
Eqs.  (7a)  and  (9),  it  was  assumed  that  T  =  T  and  u  =  u. 

The  problem  now  is  to  solve  Eqs.  (8)  and  (10a)  simultaneously  to  find  6  as 
a  function  of  The  value  of  £  at  6  =  1  will  then  define  !m  and  |3.  First, 
will  be  eliminated.  Since  6  is  a  function  of  £,  a  plot  of  6  versus  £  can 
be  made,  as  indicated  by  the  solid  line  in  Fig.  5.  The  dashed  line  in  Fig.  5 
represents  the  variation  of  6  with  |  for  a  boundary  layer  growing  under  a 
uniform  free  stream  corresponding  to  the  free  stream  condition  at  some 
station  £.  The  origin  of  this  boundary  layer,  4^,,  is  such  that  6  has  the 
correct  value  at  £.  The  boundary  layer  is  assumed  to  grow  at  £  at  the  same 
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-  ACTUAL  VARIATION  OF  8  WITH  £ . 

- VARIATION  OF  8  WITH  £  FOR  BOUNDARY  LAYER  WITH 

ORIGIN  AT  AND  FREE  STREAM  H#. 

- VARIATION  OF  8  WITH  £  FOR  BOUNDARY  LAYER  WITH 

ORIGIN  AT  £•  AND  FREE  STREAM  H#  +  AH#  . 
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rate  as  the  corresponding  uniform  free  stream  boundary  layer  (dashed  line). 
The  value  of  T  at  some  point  £  +  a£  is  then 

6  +  A6  =  JT%  +  IT)  -  £.  He  (11) 

However,  the  boundary  layer  at  £  +  A|  may  be  viewed  as  corresponding  to  a 

boundary  layer  originating  at  £.  +  A|^  and  developing  under  a  uniform  free 

stream  corresponding  to  the  free  stream  at  £  +  A£.  The  latter  free  stream 

is  characterized  by  the  value  H  +  AH  .  The  resulting  boundary  layer  is 

e  e 

denoted  by  the  dash-dot  curve  in  Fig.  5.  The  choice  of  A£.  is  such  that  the 
excess  mass  flow  in  the  boundary  layer  6  +  A^  at  £  +  A£  is  the  same  as  for 
the  boundary  layer  originating  at  Thus 


6  +  A6  =  +  Ai)  -  (4.  +  A£.)  (He  +  AH).  (12) 

Equating  Eqs.  (11)  and  (12)  and  neglecting  higher  order  terms  gives 


1  AHe 

2  4  -  il  '  He 


(13) 


which  relates  a£^  and  AHe« 

The  total  differential  of  Eq.  (10a)  is 

M  i 

6  '  2  S  -  «i  He 

In  view  of  Eq.  (13)  this  becomes 


(14) 


A6  _  1  Aj 

7  "  2  6  -  4i 


(15) 
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2  2 

But,  from  Eq.  (10a),  £-£.=£  /H  .  Hence,  Eq.  (15)  becomes 
—  —  2  l  e 

6A6/H  =  a£/2,  or  in  integral  form 


6  = 


(16) 


This  expression  can  be  integrated  to  find 
of  £  at  6 
and  (10b), 


1  is  denoted  £  and  defines 


£  as  a  function  of  5.  The  value 
I  -  In  particular,  from  Eqs.  (4) 


l>Vt*-Z‘Km*l>z  or  <17) 

The  subscript  2  is  used  when  is  obtained  from  a  numerical  integration 
of  Eq.  (16).  The  subscript  1  is  used  when  is  obtained  from  an  approxi¬ 
mate  integration,  as  will  be  discussed  later. 

In  order  to  integrate  Eq.  (16)  it  is  necessary  to  express  He  as  a  function  of 
6,  A  procedure  for  doing  this  is  outlined  by  Eqs.  (A-13)  through  (A-l6)  in 
the  Appendix.  These  equations  apply  for  an  ideal  gas  where  or  =  1  and 
pp  =  constant  across  the  boundary  layer.  A  numerical  integration  of  Eq.  (16) 
is  required.  This  has  been  done  for  y  =  7/5,  y  =  5/3,  and  values  of  W  from 
1.  100  through  6.  The  results  for  ^  are  included  in  Table  1.  These  will  be 
discussed  later. 

Approximate  Analytical  Integration  of  Eq.  (16) 

Equation  (16)  requires  a  numerical  integration  and  only  the  case  of  an  ideal 
gas  where  o*  =  1  and  pp  =  constant  has  been  evaluated.  In  the  present 
section,  an  approximate  analytical  integration  of  Eq.  (16)  is  obtained  which 
is  suggested  by  the  solution  for  strong  shocks.  It  will  be  shown  that  this 
approximate  integration  gives  accurate  results  for  all  M  ,  except  M  near  1. 

8  S 
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2 

For  M  «  1  (i.e.  ,  for  strong  shocks),  the  state  properties  of  the  fluid  in 
e  i  o 

the  free  stream  between  £  =  0  and  £  =  £  can  be  assumed  to  remain  constant. 

m  2 

That  is,  [p  /p  ]^  =  P  /p  1  +  0(M  )  .  Equations  (8)  and  (lOd)  become 

e  e,  o  e  e,  o  e,o 

6  =  1  -  V  (18a) 


He  ^  t  L)  (18b) 

I 

We  now  assume  that  1  +  I  /(f  •  *1)^  remains  fairly  constant  for  a  given  flow 
and  may  be  taken  equal  to  its  value  at  £  =  0.  (Note  from  Tables  2  and  A- 1 
that  I  /(f  -  -n)  is  small  for  strong  shocks.  However,  I  /(f  -  n)  decreases 
with  V  (Eq.  A- 1 2)  so  that  the  present  approximation  tends  to  overestimate  G.) 
Equation  (18b)  now  becomes 


H 

e 


=  n/v 


(f  -  ja 

[(f  -  tl)00]o 


G 


o 


(19a) 


W- V  «>/l  +  1. 022W  G 
W  -  1  n/v  +  1. 022W  ° 


(19b) 


where  the  substitution  for  \fv  (f  -  t|)  /( f  -  q)  in  going  from  Eq.  (19a)  to 

00  00 i o 

Eq,  (19b)  was  found  from  Eqs.  (A-8a)  and  (A- 10a).  In  effect,  Eq.  (19b) 
assumes  that  the  change  in  is  due  primarily  to  changes  in  boundary  layer 
velocity  profile.  Substitution  of  Eqs.  (18a)  and  (19b)  into  Eq.  ( 1 6)  yields 


t  = 


(1  -  V)(l  +  DV)dV 
(1  -  EV)2 


(20a) 
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where 


B  -  w  1  +  1. 022  W  2 
2.  044  _  j^2  o 

D  =  1/1. 022  W 
E  -  1  /W 

Equation  (20a)  can  be  integrated  and  gives,  for  V  =  0, 


1  I  DE  -  E  -  2D 


m  BE2 


ln(l  -  E)  -  E  -  2D 


+  +o(i/w2)] 


These  equations  become,  in  terms  of  W  and  Gq, 


(20b) 

(20c) 

(20d) 


(21a) 


(21b) 


_  6.  044  (W  -  1)' 


m 


1  +  1. 022  W 


f(w  -  rkzMirh-)  - l] 


(22a) 


-  2[i  +0(1/W2)] 


(22b) 


In  deriving  Eq.  (22b),  a  term  0.  01  5/ W  was  neglected  compared  with  1. 

Equations  (22)  can  be  used,  with  Eq.  (17),  to  find  p.  The  resulting  values 
of  p  are  given  the  subscript  1  to  indicate  that  they  are  obtained  from  the 
approximate  integration  of  Eq.  ( 1 6) .  Values  of  p^  obtained  from  Eqs.  (22), 
are  included  in  Tables  1  and  2  and  will  also  be  discussed  later. 
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The  parameter  (3^  has  a  simple  relation  to  when  W2  »  1.  Thus,  from 
Eqs.  (5),  (17),  and  (22b) 


pi  =  po-vTTl  [1  +  °(1/w2)l 


(23) 


It  is  seen  that  p^  is  larger  than  p^  but  approaches  Pq  as  W  oo. 
Simplified  Expressions  for  p^ 

Equations  (17)  and  (22)  indicate  that  p^  depends  only  on  W  and  Gq.  The 
value  of  W  is  found  from  normal  shock  relations  and  can  be  considered 
known.  The  problem  then  is  to  find  G^.  In  References  12  and  13,  analytic 
interpolation  formulas  for  Gq  have  been  presented.  These  are  used  herein 
to  obtain  simple  analytic  expressions  for  p^. 

For  an  ideal  gas  with  pp  constant  in  the  boundary  layer 


G 

o 


1.  135(W  -  1)  fj  f  ]co 

+  1. 022  W  I  (f  "  ^oJo 


(24) 


where  the  coefficient  of  the  bracketed  term  is  (f  -  r\) ^  (see  Eq.  A-8a). 
Substituting  into  Eqs,  (17)  and  (22b)  and  neglecting  a  term  0.  01 /W  compared 
with  1  yields 


P,  =  1- 59  [  1  +  1^/(1  -  n)m  ]Q  [  1  +0(1/W2)]  (25) 


The  ratio  [  l^/(i  -  ij)  1  has  been  evaluated  for  or  =  1  in  the  Appendix. 
Substituting  Eq.  (A-9b)  into  Eq.  (25)  gives 


pi  =  1,59  1  +wrf  l1  +^r~)][1  +  0(1/w2)i  <26a) 
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This  agrees  with  the  values  of  in  Table  1  for  <r  =  1  and  y  =  7/5  and  5/3 
to  within  3  percent  for  W  >  2.  A  similar  expression  can  be  derived  to 
correlate  the  cr  =  0.  72  data  in  Table  1.  If  the  constants  in  the  expression 
for  [^/(f  -  t^)  ]  are  adjusted  so  that  p^  has  the  correct  value  at  W  =  6, 

y  =  7/5  and  W  =  4,  y  =  5/3,  the  following  expression  is  obtained  for  or  =  0.  72 

This  agrees  with  the  results  in  Table  1  to  within  2  percent  for  W  >  2. 

The  effect  of  variable  pp  in  the  boundary  layer  has  been  treated  in  Reference  13 

for  the  case  of  air.  It  was  found  that  G  is  related  to  the  value  of  G  found 

o  o 

from  a  constant  pp  solution  by 

G  =  (C  )°*  37  (G  )  .  *  (27) 

o  '  e,  o'  o  pp  =  constant 

where  C  =  (p  p  /p  p  )  .  The  exponent  0.37  correlated  the  numerical 
e  i  o  e  e  w  w  o 

boundary  layer  solution  for  4  £MS<  14,  p^  =  0.001,  0.01  atmospheres,  and 
T  =  T  =  522°R  (in  Ref.  13)  and  should  be  valid  for  air  in  a  low  pressure 
shock  tube.  Equation  (27)  should  also  give  reliable  estimates  for  the  effect  of 
variable  pp  for  gases  other  than  air.t  This  leads  to  the  following  expression 

P,  =  (C  )0,  37  (p.)  .  *  (27) 

*1  e,o'  xrrpp= constant 


+  0.  37 

‘For  weak  shocks,  Ce  0  is  nearly  one  and  the  correction  for  variable  pp 
is  small.  For  strong  shocks  GQ  depends  primarily  on  (f  -  rj)^  which  is 
obtained  from  an  integration  of  the  momentum  equation  (Gf ,f) *  +  *ffH  =  0 
where  C  =  pp/pwpw.  The  quantity  C  varies  monatomically  from  1  at  the 
wall  to  CGf  o  at  the  edge  of  the  boundary  layer.  Its  effect  should  depend 
primarily  on  the  value  of  Ce>  0  and  should  not  be  a  strong  function  of  the 
nature  of  the  gas. 
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for  obtaining  p^  from  values  of  p^  found  from  constant  p \x  boundary  layer 
solutions.  Equation  (26b)  then  gives,  for  <r  =  0.  72  and  W  >  2, 


Pi 


1.  59 


c0-37 

e,  o 


\  .  0. 802 W  / 

+  2.24\] 

+  zw  -  1  V 

1  +  w  ;] 

(28a) 


where  an  ideal  gas  (constant  y)  solution  for  the  shock  is  assumed  (due  to 
the  presence  of  Z).  For  strong  shocks,  where  y  no  longer  equals  the  ideal 
value,  Z  approximately  equals  W.  Hence,  for  strong  shocks,  Z  can  be 
replaced  by  W  and  Eq.  (28a)  becomes 


PL  =  1.  59  C 


0.  37 
e ,  o 


\  A  0.  802  W  / 

2.24\ 

1  +  2  V 

-  1  V 

W 

(28b) 


Results  from  Eq.  (28b)  for  air  are  listed  in  Table  2.  These  agree  within 
2  percent  with  the  exact  values  of  p^  (computed  from  Eqs;  17  and  22a) 
for  Mg  >  4.  The  agreement  is  within  1  percent  for  Mg  >  8. 

C.  SUMMARY  AND  DISCUSSION  OF  RESULTS 

Various  estimates  for  p  have  been  made.  These  will  now  be  reviewed  and 
compared. 

The  parameter  pQ  was  found  by  assuming  the  free  stream  to  be  uniform  and 
by  finding  I  such  that  the  excess  mass  flow  through  the  boundary  layer 
equalled  the  mass  flow  through  the  shock.  This  assumes  that  the  relative 
velocity  between  wall  and  free  stream  is  u  -  u  .  Near  the  contact 
surface,  the  relative  velocity  is  actually  u  .  Hence,  pQ  is  too  small 
(i.e.  ,  I  m  is  overestimated).  For  weak  shocks  u^  is  nearly  equal  to  ug  q 
and  the  error  is  very  large. 

In  order  to  obtain  an  improved  estimate  for  p,  the  variation  of  free  stream 
conditions  was  taken  into  account  and  a  local  similarity  boundary  layer 
solution  was  employed.  A  numerical  integration  was  required  and  the 
resulting  values  of  P  were  denoted  P^.  Values  of  p^  were  obtained  for  an 
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ideal  gas  with  y  =  7/5  and  5/3,  <r  =  l,and  pp  =  a  constant  across  the  boundary 
layer.  The  results  are  given  in  Table  1  and  Fig.  6.  As  expected,  p^  is 
larger  than  Pq.  The  difference  is  very  marked  near  W  =  1  since  Pq  behaves 
like  hi W  -  1  while  P2  behaves  like  1/n/W  ”  as  W 1.  The  latter  behavior 
for  P2  is  required  in  order  that  f  be  finite  as  W  1  (see  Eq.  3).  The 
infinite  value  of  p2»  at  W  =  1,  could  have  been  avoided  if  Roshko  had  taken 
the  characteristic  velocity  to  be  u  rather  than  u  -  u  in  Eqs.  (lb)  and 
(lc). 

The  numerical  integration  to  obtain  P2  is  tedious  and  simplifications  were 
introduced  to  permit  a  closed  form  integration.  The  resulting  values  of  P 
were  denoted  p^.  Values  of  p^,  for  <r  =  1,  y  =  5/3. and  7/5,  and  pp  constant 
may  be  compared  with  the  corresponding  values  of  P2  in  Table  1  and  Fig.  6. 

It  is  seen  that  p^  is  somewhat  larger  than  P2  for  larger  W  and  is  smaller 
than  p^  for  W  near  1.  The  two  agree  within  8  percent  for  y  =  5/3  and 
5  percent  for  y  =  7/5,  except  near  W  =  1,  where  the  local  similarity  solution 
itself  is  least  accurate.  Since  p^  is  relatively  simple  to  obtain,  it  will  be 
used  henceforth  instead  of  P2  to  evaluate  the  flow  in  a  low  density  shock 
tube. 

Values  of  p1  for  air,  including  real  gas  effects  and  variable  pp,  are  given 

in  Table  2  and  Fig.  7.  The  effect  of  variable  pp  is  to  decrease  p.  by  a 

0  37  A 

factor  (C  )  This  effect  is  most  pronounced  for  strong  shocks  and 

causes  p^  to  decrease  continuously  with  Mg.  Figure  7  includes  values  of 

Pj^  (Eq.  5b),  including  real  gas  effects,  as  computed  in  References  8  and  14. 

The  agreement  between  PD  (from  Ref.  8)  and  P,  near  M  =  10  is  fortuitous, 

xv  IS 

being  due  to  a  rough  approximation  by  Roshko  for  the  effect  of  C  .  The 

e ,  o 

results  of  Reference  14,  for  p^,  are  based  on  an  accurate  boundary  layer 

solution.  It  is  seen  that  p^  underestimates  p^  and  therefore  considerably 

overestimates  t  [recall  I  ~  P  ]  for  the  range  of  M  in  Fig.  7. 
m  m  s 

Values  of  p.  for  argon  are  included  in  Fig.  7.  These  were  obtained  by 

1  0  37 

multiplying  the  values  of  p.  in  Table  1  (or  =  0.  72,  y  =  5/3)  by  (C  ) 

i  e ,  o 

The  normal  shock  solution  for  argon  was  based  on  y  =  5/3  and  the  results 
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.000  oo  --  1.60  1.89  --  --  1.53  1.81  1.73 


T  Values  of  P  and  P .  are  obtained  from  data  in  Table 


Fig.  7.  Values  of  p  for  Air  (Based  on  Equilibrium  Real  Gas  Properties, 

p  »  0.  001  atm)  and  for  Argon,  Including  Effect  of  Variable  pp; 

°°  T  =  522°R 

oo 
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are  therefore  not  extended  beyond  Mg  =  10  (where  ionization  effects  become 
important) . 

0 

Roshko  measured  test  time  in  a  low  density  shock  tube  for  2.  5  <  Mg  <  9  in 
air  and  argon.  He  found  that  p  =  n/3  gave  a  mean  correlation  of  the  experi¬ 
mental  data.  However,  since  he  computed  p  from  Eqs.  (5),  he  had  no 
rational  explanation  of  why  p  should  be  so  high.  The  present  results  (Fig.  7) 
indicate  that  p  -  is  a  reasonable  mean  value  for  his  tests.  The  marked 
decrease  of  p  with  Mg  indicates  that  smaller  Mach  number  ranges  should 
be  used  when  correlating  the  experimental  data,  particularly  for  moderate 

M  . 
s 

The  present  results  for  p.  have  been  used  in  Eq.  (3)  to  obtain  the  variation 
2  1 

of  I  p  ./d  p  with  M  for  air  and  argon.  The  results  are  given  in  Fig.  8. 
m  st  *  oo  s 

The  standard  condition  in  Eq.  (3)  was  taken  to  be  Tg^  =  522  R  and 
pgt  =  1  atmosphere,  for  which 


(y).t  =  6-93xlo6ft‘1 

air 

=  7.  39  X  106  ft"1 

argon 

Roshko  used  the  values  6.  74  X  10^  and  7.  03  X  10^  for  air  and  argon, 

respectively.  These  values  correspond  to  a  higher  standard  temperature 

[Note  that  (pa/p)^  *  Tgt”*‘  ^  for  air,  since  p  ~  T^*  The  effect  of 

initial  temperatures  other  than  522°R  can  be  readily  taken  into  account]. 

The  results  for  air  in  Fig.  8  are  based  on  values  of  p^  obtained  from 

Fig.  7  for  M  <  14  and  on  values  obtained  from  Eq.  (28b).  Real  gas 

normal  shock  solutions  were  employed.  The  difference  in  the  values  of 

I  for  p  =0.5  cm  and  p  =0.  001  cm  is  due  to  the  effect  of  initial 
m  rco 

pressure  on  the  shock  solution.  This  difference  is  small  and  probably  within 
the  accuracy  of  the  present  solution.  The  results  for  argon  in  Fig.  8  are 
based  on  the  values  of  Pj  in  Fig.  7  and  an  ideal  gas  normal  shock  solution 
(V  =  5/3). 
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III.  UNSTEADY  SOLUTION 


In  the  previous  section,  the  asymptotic  value  of  the  separation  distance 
between  the  shock  and  contact  surface  was  found.  This  defines  the  test  time 
in  shock  tubes  that  are  sufficiently  long  to  permit  the  asymptotic  steady  state 
flow  to  be  reached.  The  separation  distance  as  a  function  of  time  is  required 
in  order  to  estimate  test  times  in  shock  tubes  wherein  the  steady  state  solution 
is  not  achieved.  This  problem  is  now  discussed. 

It  is  desired  to  find  the  separation  distance  as  a  function  of  time,  t,  assuming 
that  the  shock  and  the  contact  surface  are  coincident  at  t  =  0.  This  problem 
is  very  difficult  to  solve,  and  it  is  necessary  to  set  up  an  approximate  model 
for  the  flow.  Two  limiting  approaches  seem  tractable.  The  first  is  to  assume 
that  the  contact  surface  moves  with  uniform  velocity  and  to  attempt  to  compute 
the  nonuniform  motion  of  the  shock  relative  to  the  contact  surface.  At  first 
glance,  References  8  and  10  appear  to  use  this  approach.  'However,  in  both 
of  these  references  it  is  assumed  that  the  mass  flow  through  the  shock,  rr^, 
is  a  constant  and  thus  a  uniform  shock  velocity  is  implicitly  assumed.  In 
addition,  the  proper  boundary  layer  theory  to  be  used  in  the  contact  surface 
fixed  co-ordinate  system  is  not  clear.  (In  Refs.  8  and  10  attempts  are  made 
to  relate  the  boundary  layer  to  that  behind  a  constant  velocity  shock  with  a 
uniform  free  stream  downstream  of  the  shock.  )  The  second  approach  is  to 
assume  that  the  shock  moves  with  uniform  velocity  and  to  find  the  rate  at 
which  the  contact  surface  moves  relative  to  the  shock.  This  model  can  be 
specified  more  precisely  than  the  first  one,  is  self-consistent,  and  leads 
naturally  to  the  asymptotic  solution  obtained  in  the  previous  section.  The 
latter  model  will  be  treated  herein. 

A-  SEPARATION  BETWEEN  SHOCK  AND  CONTACT  SURFACE 

Consider  flow  in  a  co-ordinate  system  in  which  the  shock  is  fixed  and  the 
wall  moves  with  velocity  u  .  Assume  that  at  time  t  =  0  the  contact  surface 
coincides  with  the  shock  and  that  at  some  later  time,  t,  the  portion  of  the 
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contact  surface  which  is  in  the  free  stream  is  located  at  I  (Fig.  9).  Also 
assume  that  the  flow  between  the  shock  and  the  contact  surface  is  steady. 

For  steady  flow,  the  rate  of  mass  flow  through  the  shock  equals  the  rate  of 
mass  flow  through  a  control  surface  at  t  .  Thus 

(p  u  )  A  =  fp u  dA 
re  e  o 

=  p  u  A  +  p  u  L(-6  )  (29) 

re  e  re  e 

jje 

where  6  is  the  displacement  thickness  at  t .  This  continuity  equation  can 
be  reduced  to 


S  = 


peUe 

(p  u  ) 
e  e  o 


(30) 


which  is  the  same  as  Eq.  (8).  If  f  is  replaced  by  £  (Eq.  10b),  and  local 
similarity  is  assumed  for  the  boundary  layer,  it  again  follows  (see  derivation 
of  Eq.  16)  that 


d£  =  2Td7/H2 

e 


(31) 


We  are  now  interested  in  obtaining  I  as  a  function  of  t.  But  df  /dt  =  u  . 

<  e 

Thus  t  =  /  dl/u  .  Define  X  =  u  t/f  ,  T  =  I /l  It  then  follows 
o  e  e ,  o  m  m 

(from  t/f  =  i/i  and  Eq.  31)  that 
mm 


X 


u  t 
e,  o 


m 


x 

s 


/o^dlT/IVH*) 

~PTd¥/Hr~ 

Jo  e 


(32a) 
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Fig.  9.  Flow  Model  for  Finding  Separation  Distance  Between 
Shock  and  Contact  Surface  as  a  Function  of  Time 
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The  quantity  xg  in  Eq.  (32a)  is  the  distance  of  the  shock  from  the  diaphragm, 
xg  =  u^t.  Equations  (32)  give  the  separation  distance  !  as  a  function  of 
t  (or  xg).  The  parameters  X  and  T  are  the  same  as  those  introduced  by  Roshko, 
The  quantity  T  is  the  ratio  of  the  separation  distance  f  to  the  asymptotic 
value  |  .  The  quantity  X  is  the  time  t  divided  by  the  time  it  would  take  an 

ideal  shock  tube  (no  boundary  layer)  to  achieve  a  separation  distance  |  . 

The  latter  may  also  be  viewed  as  the  shock  distance  xg  divided  by  the  distance 
required,  in  an  ideal  shock  tube,  to  achieve  a  separation  distance  fm>  Note 
also  that  dT/dX  =  u^/u^  q.  For  an  inviscid  flow,  X  =  T,  The  departure  of  X 
from  T  is  a  measure  of  the  wall  boundary  layer  effect. 

Equations  (32)  have  been  integrated  for  an  ideal  gas  with  cr  =  1  and  constant 
pp  across  the  boundary  layer  for  y  =  7/5  and  5/3.  Equations  ( A-  1  3)  through 
(A- 16)  were  used.  Plots  of  X  versus  T  are  given  in  Fig.  10.  The  results 
for  W  =  1.25,  y  =  7/5,  and  for  W  =  1.25,  y  =  5/3,  are  sufficiently  close  that 
they  appear  as  a  single  curve  in  Fig.  10.  Similarly,  the  results  for  W  =  6, 
y  =  7/5,  and  W  =  4,  y  *  5/3,  appear  as  a  single  curve.  The  displacement 
between  these  curves  is  small  and  intermediate  values  of  W  are  therefore 
not  plotted. 


A  simpler  gelation  between  X  and  T,  suggested  by  the  solution  for  strong 
shocks,  can  be  found  as  follows.  Write  Eq.  (29)  as 


1 


u  ) 
e  o 


Pe  ue  L(-6#) 

~Tp  "  u  J  “A 
e  e  o 


(33) 
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Fig.  10.  Variation  of  T  With  X  From  Numerical  Integration  of  Eq.  (32) 
(For  Ideal  Gas,  <r  =  1,  and  pp  =  Constant)  and  From  Eq.  (35) 


If  variations  in  free  stream  density  are  neglected,  the  first  term  on  the  right 
hand  side  can  be  written  as  ^e/^e  Q  =  dT/dX.  The  second  term  on  the  right 
hand  side  is  the  ratio  of  the  excess  mass  flow  through  the  boundary  layer  at 
X  to  the  mass  flow  through  the  boundary  layer  at  f  .  For  strong  shocks, 
the  external  free  stream  remains  fairly  uniform  and  the  excess  mass  flow  in 
the  boundary  layer  grows  approximately  as  *JT  so  that  the  second  term  can  be 
approximated  by  n/i f7F^  =  n/T\  Equation  (33)  then  becomes 

1  =  +  n/t  (34) 

This  is  the  same  equation  as  that  derived  by  Roshko  and  the  solution  is 

-j  =  t  n  (1  -  sff)  +  n/t  (35) 

Equation  (35)  is  also  plotted  in  Fig.  10  and  is  indistinguishable  from  the 
previous  results  for  W  =  6,  y  =  7/5,  and  W  =  4,  y  =  5/3.  This  shows  that 
the  plot  of  X  versus  T  is  relatively  insensitive  to  W  and  y  and  that  Eq.  (35) 
gives  an  accurate  representation  except  for  W  very  near  1. 

B.  TEST  TIME 

In  the  previous  calculation,  the  separation  distance  was  obtained  as  a  function 
of  t.  A  quantity  that  is  perhaps  oi  greater  interest  is  the  test  time  (i.  e.  ,  the 
difference  in  time  between  the  arrival  of  the  shock  and  the  arrival  of  the 
contact  surface)  at  a  fixed  value  of  x.  This  quantity  will  now  be  discussed. 

Designate  the  test  time  by  t.  For  t-*oo,  t  =  {  /u  .  Define  T  =  u  r/l  4 

7  m  w  w  , 

which  is  the  test  time  at  x  divided  by  the  test  time  at  x-*oo.  Since  u  -  u 

w  e 

is  the  velocity  of  the  contact  surface  relative  to  the  shock  tube  wall,  the  test 

time  will  nearly  equal  t  /  (u  -u  ).  Also,  u  >  u  >0.  It  then  follows  that 

we  e  ,o  —  e  — 

W  — 

flTT-i  T  >  T  >  T  (36) 
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At  the  start  of  the  motion,  t  =  WT/(W  -  1),  whereas  after  long  times  r  =  T. 

If  W  is  large,  t  =  T  throughout  the  entire  motion. 

The  test  time  t  at  a  particular  station  x  is  found  analytically  as  follows, 
a  a 

In  the  present  model,  the  shock  moves  with  uniform  velocity  u^,  Let  I  be 
the  separation  distance  when  the  shock  is  at  x^,  and  let  x^  be  the  location  of 
the  shock  when  the  contact  surface  is  at  x  ,  From  Fig.  11  it  is  clear  that 

a 

x  =  x-  -  I-  and  t  =  t,  -  t  =  I  ,/u  .  In  nondimensional  variables,  the 
a  d  b  a  b  a  b  w 

latter  expressions  become 

x,  *  *b  -  Tb/w  (37*» 


Equations  (37)  give  corresponding  values  of  X  and  t  from  corresponding 
values  of  X^  and  (which  are  found  from  Eqs.  (32)  or  Fig.  10).  Since 
station  a  is  arbitrary,  the  subscript  a  may  be  removed  from  Eqs.  (37).  Plots 
of  t  versus  X  are  given  in  Fig.  12  for  an  ideal  gas,  where  <r  =  1  and 
p  p  =  constant. 

If  Eq.  (35)  is  used,  the  relation  between  X  and  T  becomes 

-X=  2  [in  (1  -  s/t)  +  n/T]  +T/W  (38) 

Equation  (38)  is  also  plotted  in  Fig.  12. 

It  is  seen  from  Fig.  12  that  the  variation  of  t  with  X  is  more  sensitive  to  W 
than  is  the  variation  of  T  with  X  (Fig.  10.  )  Equation  (38)  gives  a  reliable 
estimate  for  the  variation  of  t  with  X  except  for  W  near  1. 
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Fig.  12.  Test  Time  at  Given  Station 


') 

Since  the  variation  of  T  with  X  is  relatively  insensitive  to  W  (compared  with 
t  versus  X),  it  is  preferable  to  attempt  to  correlate  experimental  data  on  the 
basis  of  T  versus  X.  The  data  should  then  correlate  with  the  single  curve 
defined  by  Eq.  (35)  for  all  values  of  W  except  W  near  1.  Experimental  obser¬ 
vations  of  test  time  can  be  reduced  to  T  as  follows.  If  t  is  the  nondimension- 

a 

alized  test  time  at  X  ,  Eq.  (37)  indicates  that  the  corresponding  values  of  T 
a 

and  X  are 


X  =  X  +  T  /W 

(39a) 

a  a 

T  =  T 

a 

(39b) 

This  reduction  procedure  becomes  more  important  as  W  decreases.  A 
similar  procedure  has  been  discussed  in  Reference  10. 
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IV.  COMPARISON  WITH  REFERENCES  8  AND  10 


References  8  and  10  used  a  contact  surface  fixed  co-ordinate  system  to  study 
the  separation  between  the  shock  and  the  contact  surface.  The  derivation  of 
their  basic  equations  is  outlined  here.  The  assumptions  which  are  required 
to  reduce  the  equations  for  the  shock  fixed  model  to  those  of  References  8 
and  10  are  also  noted  and  discussed. 

The  notation  to  be  used  in  the  contact  fixed  co-ordinate  system  is  indicated 
in  Fig.  13.  Subscripts  1  and  2  denote  conditions  upstream  and  downstream 
of  the  shock,  respectively.  Ufi  is  the  velocity  of  the  shock  relative  to  the 
wall  and  u^  is  the  velocity  of  the  free  stream  in  region  2  relative  to  the  wall. 

In  both  References  8  and  10  it  is  assumed  that  the  free  stream  is  uniform  in 
region  2  so  that  u^  is  a  constant  and  also  equals  the  velocity  of  the  contact 
surface  relative  to  the  wall.  Hence,  the  free  stream  velocity  is  zero  in 
region  2  in  contact  surface  fixed  co-ordinates.  The  above* notation  is  the 
same  as  that  used  in  References  8  and  10.  In  addition,  the  velocity  in  the 
boundary  layer,  relative  to  the  contact  surface,  will  be  denoted  by  &  (Fig.  13). 

The  mass  balance  between  flow  entering  and  leaving  region  2  is  found  as 
follows.  The  rate  at  which  mass  enters  the  shock  at  any  instant  is 

=  P1  Us  A  =  p2  <Us  ‘  u2*  A  <40) 

The  rate  at  which  mass  leaves,  at  the  plane  of  the  contact  surface,  is 

mc  =  Lj£  pu  dy  (41) 

where  pu  is  zero  outside  of  the  boundary  layer.  The  rate  of  increase  in  mass 
between  the  shock  and  contact  surface,  rh,  can  be  found  in  the  following 
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(a)  TIME  t 


(b)  TIME  t  +  dt 

Fig.  13,  Flow  in  Contact  Surface  Fixed  Coordinate 
System  of  References  8  and  10 


manner.  Consider  the  interval  of  time  between  t  and  t  +  dt.  The  separation 
distance  increases  by  an  amount  df  .  Assume  that  the  mass  between  0  and  f 
at  time  t  is  the  same  as  the  mass  between  df  and  I  +  d|  at  time  t  +  dt  (see 
Fig.  13).  The  increased  mass  between  the  shock  and  contact  surface  is  then 
contained  between  stations  0  and  df  and  equals 

dm  =  (fp  dA)  di 

=  [  p2A  +  L  fQ  (p  -  P2)  dy]  df  (42) 

This  g  ives  the  following  expression  for  the  rate  of  increase  of  mass  in 
region  2: 

rii  =  [  p2  A  +  L^°°(p  -  p2)  dy]  (43) 

(Roshko  used  rh  =  A  df  /dt  and  Hooker  used  the  integral  in  the  form 

p  dy.  Reference  14  corrected  Hooker’s  integral.  )  Continuity  of  mass, 

-  rhc  =  rh,  then  gives 

IU»  -  u2)  -  >/T  (j  Jo  |j^f  )  =  Ro.hko  (44.) 

*  af  [‘  *  '  ‘Xs1)]  Hook"  (44b) 

where  the  integrations  are  te(ken  at  the  contact  surface. 
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Since  the  free  stream  was  assumed  uniform  in  region  2,  the  integrals  in 
Eq.  (44)  can  be  evaluated  using  the  results  of  References  11  to  13  and  are 
independent  of  f  and  t.  Hence  Eqs.  (44)  can  be  integrated  to  find  I  as  a 
function  of  t.  In  integrating  Eq.  (44),  both  Roshko  and  Hooker  assumed 
that  U#  -  U-,  is  a  constant.  However,  from  Fig.  13,  it  is  seen  that 
Ug  -  u2  =  <*l  /dt  and  therefore  cannot  be  a  constant  for  shock  tube  flows 
when  the  wall  boundary  layer  effects  are  important.  Hence  there  is  an 
inconsistency  in  this  model  beyond  the  assumption  of  uniform  flow  in  region  2. 
The  assumption  that  U  -  u,  is  constant  is  equivalent  to  assuming  that  the 
mass  flow  rate  through  the  shock  is  constant.  This  is  precisely  the  basic 
assumption  made  in  the  shock  fixed  model  and  it  is  not  unexpected  that  the 
two  models  give  a  similar  variation  of  I  with  t. 

Equations  (44)  can  be  put  into  shock  fixed  notation  by  letting  U  -  u,  =  u  , 

s  &  e ,  o 

vl  =  u  -  u  ,  and  p7  =  p  The  result  is 

e ,  o  c»  e ,  o 


e,  o 


d! 

HF 


00 


(— 

\v, 


e,  o'  e,  o 


Roshko  (45a) 


u 

e, 


o 


df  L 
ST  X 


e,  o 


Hooker  (45b) 


(Roshko  used  pu/p^  =  pu/pe  Q  -  ue  q  as  well  as  pu/ =  p(u  -  u^  Q)/Pe  0 
and  the  former  is  the  transformation  error  previously  noted  in  connection 
with  Eq.  5b.  )  The  limit  di/dt  =  0,  in  Eqs.  (45),  shows  that  the 
asymptotic  separation  distance  I  is  obtained  from  a  value  of  (3 
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equal  to  ae  defined  in  Eq.  (5b).  When  put  in  term*  of  X  and  T,  Eq.  (45a) 
reduces  to  Eq.  (34)  and  the  integral  ie  given  by  Eq.  (35).  However,  in  terms 
of  X  and  T,  Eq.  (45b)  becomes 


.-8"* 


■W. 


*8 


m 


which  can  be  integrated  to  give 


•7  *  <47> 

Equation  (47)  is  Hooker’s  modification  of  Roshko's  Eq.  (35).  The  difference 
between  these  equations  is  small,  particularly  for  strong  shocks. 

Equations  (45)  can  be  compared  with  the  corresponding  equations  which 
arise  in  the  shock  fixed  model.  Equation  (29)  can  be  written 

1  ‘  Wo  ‘  **  Wo  Ct  (t  '  ‘)# 


p  u 
re  e 


Wo^*/ 


Of- 


(48) 


This  can  be  reduced  to  Eqs.  (45)  by  taking  p  =  p  and  by  taking  u  equal 

to  either  u  or  df  /dt,  depending  on  the  particular  term  involved.  If 
e ,  o 

u  =  u  in  the  third  term  of  Eq.  (48)f  and  if  u  =  df  /dt  in  the  forth  term, 
c  e  ,  o  e 

Eq.  (48)  reduces  to  Hooker's  Eq.  (45b).  However,  it  is  more  consistent 

to  let  u  =  u  for  both  these  terms.  This  would  mean  that  the  coefficient 
e  e,  o 

(df  /dt)/u  on  the  right  hand  side  of  Hooker's  Eq.  (45b)  could  be  replaced 
e  ( o 
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by  1.  If  this  substitution  is  made  in  Eq.  (45b)  then  it  would  yield  a  value  of 
P  defined  by  Pq  (which  is  more  accurate  than  PR)  and  a  variation  of  T  with  X 
as  given  by  Roshko's  Eq.  (35).  The  latter  was  shown  in  Fig.  10  to  be  in  good 
agreement  with  the  integral  solutionof  the  shock  fixed  model  for  all  W  except 
W  near  1. 

Hence  Hooker's  relatively  small  correction  (Eq.  47)  to  Roshko's  result 
(Eq.  35)  for  the  variation  of  T  with  X  doeB  not  appear  warranted.  Equation  (35) 
is  sufficiently  accurate,  considering  the  basic  limitations  of  the  contact 
surface  fixed  model.  The  major  problem  is  that  of  accurately  determining  p. 


V.  COMPARISON  WITH  EXPERIMENTAL  DATA 


Experimental  observations  of  test  time  in  low  pressure  shock  tubes  are 
reported  in  References  7,  8,  10,  and  15.  The  results  of  the  first  three  of 
these  papers  are  summarized  in  Fig.  14.  The  results  will  be  considered  in 
the  light  of  the  present  study.  In  particular,  the  question  of  whether  the  use 
of  leads  to  improved  correlations  will  be  investigated. 

The  experimental  results  in  Fig.  14  were  reduced  on  the  basis  of  p  =  s!3 

and  T  =  t  .  ^  The  value  of  P  =  was  chosen  by  Roshko  as  a  mean  fit  to 

his  experimental  data.  Recall  that  both  X  and  T  are  multiplied  by  p^  (since 
-2 

f  ~  P  ).  An  increase  in  P  will  result  in  an  increase  in  both  X  and  T.  In 
regions  where  T  has  reached  its  asymptotic  value  (i.  e.  ,  is  independent  of  X), 
changes  in  p  will  raise  of  lower  T,  but  the  effect  on  X  will  be  unimportant 
with  regard  to  correlation  with  the  asymptotic  portion  of  the  theoretic  h) 
curve  in  Fig.  14.  If  p  were  correctly  evaluated,  the  experimental  data  should 
by  asymptotic  to  T  =  1.  Hence,  the  asymptotic  data  can  be  used  to  evaluate 
P.  In  particular,  P  =  n/JTt^  where  T^  is  the  asymptotic  value  of  T  in  Fig.  14 
for  the  data  in  question. 

The  results  of  the  above  procedure  are  compared  in  Table  3  with  values  of  P^ 
from  Fig.  7.  For  air  in  the  range  Mg  =  5  through  9,  Pj  is  too  small  by  about 
10  percent.  The  results  for  argon  cover  three  Mach  number  ranges  and  are 
more  complete.  It  is  seen  that  Pj  is  too  small  by  about  35  percent  for 
M  =1.6,  20  percent  for  M  =4,  and  10  percent  for  M  =5  through  9.  These 
figures  are  only  approximate,  due  to  uncertainty  and  scatter  in  the  experimental 


The  quantity  t  was  found  from  an  experimental  measurement  of  test  time. 

The  procedure  noted  in  Eqs.  (39)  should  be  used  to  reduce  the  data  to  obtain 
a  plot  of  T  versus  X.  This  would  displace  the  points  in  Fig.  14  to  the  right. 
This  displacement  is  unimportant  in  the  region  where  T  is  near  its  asymptotic 
value,  which  is  the  region  of  primary  concern  here. 
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data.  However,  the  trend  appears  reasonable  in  view  of  the  fact  that  the 
procedures  used  to  obtain  Pj  and  p^,  in  the  present  analysis,  become  more 
accurate  as  increases.  Mixing  and  diffusion  at  the  contact  surface  would 
also  tend  to  yield  experimental  values  of  p  which  are  somewhat  larger  than  p^. 
The  magnitude  of  p^,  and  its  variation  with  Mg,  is  in  much  better  agreement 

than  would  be  obtained  by  using  p_.  The  large  values  of  p,  as  M  -*1,  indi- 

s 

cate  that  in  the  low  range  the  test  time  is  much  less  than  previous  esti* 
mates  based  on  p„. 

x\ 

On  the  basis  of  the  above  data,  it  appears  that  the  present  estimates  for  p^ 

are  correct  to  about  10  percent  for  M  >5.  Additional  experimental  data 

s 

are  required  to  better  define  p.  A  narrow  range  of  Mg  must  be  used  when'p 
is  evaluated  experimentally  in  the  low  range. 


Table  3.  Evaluation  of  p  from  Experimental  Data  in  Fig.  14 


Gas 

Ref. 

M 

8 

ta 

(Approx. ) 

P  =>/3/Ta 

Pi  ! 

(Fig.  7) 

Argon 

7 

'  1.6 

0.  2 

3.9 

2.4 

Argon 

10 

4 

0.6 

2.2 

1.8 

Argon 

8 

5-9 

0.9 

1.8 

1.7-  1.4 

Air 

8 

5-9 

1.3 

1.  5 

1.6  -  1.2 
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VI.  ADDITIONAL  CONSIDERATIONS 


The  present  theory  can  be  used  to  estimate  the  flow  nonuniformity  in  low 
pressure  shock  tubes.  In  addition,  the  asymptotic  shock  strength  for  given 
initial  shock  tube  conditions  can  be  estimated.  Since  in  the  present  theory  a 
laminar  wall  boundary  layer  is  assumed,  it  is  also  of  interest  to  discuss  tran¬ 
sition  to  a  turbulent  boundary  layer.  Finally,  it  has  been  assumed  that  the 
boundary  layer  is  thin  compared  with  the  tube  diameter,  and  this  assumption 
requires  verification.  These  topics  will  now  be  briefly  considered. 

A.  FLOW  NONUNIFORMITY 

An  estimate  of  the  flow  nonuniformity  between  the  shock  and  the  contact  surface 
can  be  obtained  from  Eq.  (34),  which,  in  dimensional  variables,  can  be  written 


u 


u 

e 


e,  o 


m 


(49) 


This  gives  the  variation  of  ug  with  distance  behind  the  shock.  Isentropic  flow 
relations  can  then  be  used  to  find  the  variation  with  i  of  all  other  flow  prop¬ 
erties.  In  shock  fixed  co-ordinates,  the  flow  downstream  of  the  shock  cor¬ 
responds  to  subsonic  flow  in  an  expanding  channel.  The  density,  temperature, 
and  pressure  increase  with  distance  behind  the  shock.  The  net  change  in  these 
properties  can  be  expressed  in  terms  of  the  flow  Mach  number  directly  behind 

the  shock  u  /a  and  is  small  for  strong  shocks.  The  present  discussion 
e ,  o  e ,  o 

neglects  relaxation  phenomena  which  may  also  contribute  to  variations  in  the 
free  stream  properties. 

B.  COMPARISON  OF  IDEAL  AND  ACTUAL  ASYMPTOTIC  SHOCK  STRENGTH 

The  problem  here  is  to  estimate  the  actual  asymptotic  shock  strength  that  will 
result  from  given  initial  conditions  in  a  low  pressure  shock  tube.  Perfect  gases 
with  constant  specific  heats  are  considered.  Here  the  conventional  subscripts 
1  and  4  are  used  to  indicate  initial  conditions  in  the  driven  and  driver  sections, 
respectively. 
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If  wall  boundary  layer  effects  are  neglected,  the  resulting  ideal  shock  Mach 
number  .  is  related  to  the  initial  pressure  ratio  P^/Pj  and  sound  speed 
ratio  a^/a^  by 


Z.  -  1  a.  M2  .  -  1 

l  1  S,  X 

TA1a-\)——$T  7" 

r  4  '  4  s, i 


/'».  <zi +  -  m 

>  z>  ) 


i/(z4+i) 


(50  a) 


In  the  derivation  of  Eq,  (50a),  it  is  assumed  that  the  flow  is  uniform  between 
the  shock  and  contact  surface. 

Duff  has  pointed  out  that  the  ideal  shock  tube  equations  can  be  modified  to 

account  for  the  character  of  the  asymptotic  flow  in  a  low  pressure  shock  tube. 

In  particular,  Duff  neglects  the  boundary  layer  effect  in  the  driver  gas  (valid 

for  strong  shocks).  He  also  assumes  that  the  contact  surface  moves  with  the 

same  velocity  as  the  shock  and  that  the  pressure  at  the  contact  surface  equals 

the  stagnation  pressure  (relative  to  the  shock)  of  the  flow  directly  behind  the 

shock.  The  shock  Mach  number  M  A  that  then  results  from  given  initial 

7  S ,  A 

conditions  is  found  from 


M 


uAi.i/; 

^  “4  \I 

X  + 


pi  (zi  ♦  ‘)m;.a 

p4  Tx 


4 


/(Z4+l) 


Ms,A+Zl+  1 


(z i  -  m<z1+  dm‘iA 


(Z,  +  1)/2(Z4+1) 


(50b) 


The  quantity  M  A  may  be  termed  the  actual  asymptotic  shock  strength. 

Eq,  (50b)  differs  in  form  from  Eq.  (4)  in  Reference  7,  since  Duff  was  primarily 

interested  in  determining  the  initial  pressure  ratio  P^/pj  required  to  achieve 

a  given  shock  strength  M  .  * 

s , 
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When  the  shock  tube  diaphragm  first  bursts »  the  boundary  layer  effects  are 
small  and  the  shock  Mach  number  will  tend  to  be  M  ,  as  given  by  Eq.  (50a). 

8  |  l 

Hence,  M  .  may  also  be  viewed  as  the  initial  shock  Mach  number.  (The 

8,1 

finite  time  required  to  rupture  the  diaphragm  tends  to  alter  this  result  some¬ 
what).  If  the  flow  reaches  its  steady  asymptotic  limit,  the  shock  Mach  number 
is  given  by  Eq.  (50b).  These  two  Mach  numbers  can  be  readily  compared  for 
large  initial  pressure  ratios.  In  the  latter  case,  Eqs.  (50a)  and  (50b)  can  be 
equated  to  yield 


M 


TT 


s,  A 


s ,  i 


W.  -  1 

1 

i  L 


+  °(p1/p4) 


1/(Z 


,♦» 


(51) 


where  W.  is  the  density  ratio  across  the  shock  corresponding  to  M  For 
i  s ,  1 

strong  shocks,  M  A  is  only  slightly  less  than  M  The  difference  becomes 

S  ,  A  S  ,  1 

more  pronounced  for  weaker  initial  shocks. 

g 

Roshko  has  deduced  an  expression  similar  to  Eq.  (51)  from  a  simple  physical 

argument.  Namely,  for  ideal  flow,the  shock  velocity  relative  to  the  wall  is 

U  .  =  (W  u  ).  .  For  the  asymptotic  flow,  the  shock  velocity  approximately 
s ,  i  e ,  o  i 

equals  the  ideal  contact  surface  velocity  relative  to  the  wall;  U  =  (U  -  u  ).. 

s ,  A  s  e,  o  1 

The  latter  expressions  give  Eq.  (51). 

C.  TRANSITION 

It  has  been  assumed  throughout  that  the  wall  boundary  layer  is  laminar  between 
the  shock  and  the  contact  surface.  This  assumption  is  generally  valid  for  low 
pressure  shock  tubes.  It  is  of  interest  to  compute  the  Reynolds  number  at  the 
contact  surface  for  the  asymptotic  flow  condition.  If  this  is  below  the  transition 
Reynolds  number,  then  the  validity  of  the  laminar  flow  assumption  is  established. 


-53- 


In  Reference  16,  a  transition  Reynolds  number  defined  by 


(Re)t 


u  [W  -  l]2l. 
e,oL _ J  t 

v 

e ,  o 


(52) 


has  been  proposed  for  correlating  transition  in  shock  tubes.  Here,  i  is  the 
distance  between  the  shock  and  the  transition  point.  (The  characteristic  velocity 
used  to  deduce  Eq.  (52)  is  the  flow  velocity  relative  to  the  wall  and  the  charac¬ 
teristic  distance  is  the  distance  a  particle  moves  in  the  free  stream  from  the 
instant  it  is  set  into  motion  to  the  instant  it  is  at  the  transition  point.  )  Experi¬ 
mental  observations  of  (Re)  are  summarized  in  References  16  and  17.  For 

6  * 

weak  shocks,  (Re)^  =  0(10  ).  The  value  of  (Re)^  increases  with  particularly 

for  strong  shocks  where  the  low  wall-to-free-stream  temperature  ratio  tends  to 

stabilize  the  boundary  layer.  This  stabilization  seems  to  occur  at  about 

T  /T  0.  1,  which  corresponds  to  M  %  10  in  air.  The.  data  in  Reference  17 

w  e,  o  s 

indicates  that  for  1  <  M  <9,  (Re)  appears  to  be  in  the  range 

-6  8  * 

0.  5  <  (Re)  X  10  <4.  For  larger  M  the  values  of  (Re)  tend  to  increase 

t  7  s  7  t  j  7 

markedly  and  values  as  high  as  10  and  5X10  have  been  observed  for 
around  10.  An  analytical  study  of  the  stabilizing  effect  of  low  T^/T^  q  has  been 
presented  in  Reference  18.  The  latter  study  is  consistent  with  the  experimental 
data  for  very  weak  shocks,  but  indicates  infinite  stability  for  >  2.  18  in  air. 


The  Reynolds  number  at  the  contact  surface,  corresponding  to  a  shock-contact 
surface  separation  distance  of  f  ,  can  be  used  to  determine  whether  or  not  the 
laminar  boundary  layer  assumption  is  correct.  This  Reynolds  number  is 
(Re)|  ^  =  Ue  otW  ~^^rr/Ve  0  as surning  a  unifor m  free  stream ,  and  can  be  put 


m 

in  the  form 
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( 


i 


The  right  hand  side  depends  primarily  on  M  .  For  a  given  M  ,  (Re)j  varies 

2  t  8  *  m 

a«  (dPco/pst)  .  If  d  is  increased  to  increase  then  smaller  values  of 

may  be  required  in  order  to  assure  laminar  flow.  Eq.  (53)  has  been  evaluated! 

employing  the  same  data  as  used  to  obtain  Figs.  7  and  8,  and  the  results  are 

presented  in  Fig.  15. 

The  present  results  have  been  used  to  estimate  (Re)|  for  one  series  of 

m 

Roshko's  tests  in  air.  These  tests  were  made  in  a  tube  with  d  =  1/6  feet, 

5  <  M  <9  and  0.  1  mm  Hg  <  p  <5  mm  Hg.  Using  a,n  average  value 

(pgt/poo)^(Re)/  m/d^  =  2  X  10*^°{from  Fig.  15),  it  is  found  that 

10^  <  (Re)#  <  2  X  107.  At  the  higher  values  of  (Re)#  ,  transition  to  a  tur- 

m  —  z  m 

bulent  boundary  layer  might  have  occurred.  This  would  reduce  and  increase 

the  effective  value  of  p,  as  compared  with  a  completely  laminar  boundary  layer. 

D.  BOUNDARY  LAYER  THICKNESS  AT  i 

m 

Let  6  represent  the  value  of  y  at  which  (u  -  u)/(u  -  u  )  =  0.  99  and  let 

u  w  v/  e ,  o 

be  the  corresponding  value  of  r\  .  Hence,  6^  is  a  measure  of  the  boundary 
layer  thickness.  Assuming  that  the  flow  behind  the  shock  is  uniform  6^  is  given 
by  (e.  g.  ,  Ref.  13) 


The  value  of  6  at  i  can  be  used  as  an  index  to  determine  whether  the 
u  m 

boundary  layer  thickness  is  small  relative  to  the  tube  hydraulic  radius  d/2, 
Combining  Eqs.  (2)  and  (54)  gives 


(n  - 1  ) 

v  'u  OO  P 

p*/2(  W-~l) 


(55) 


tThis  dependence  is  due  to  the  fact  that  im  varies  as  d^Poo/p»t  an<*  Reynolds 
number  per  unit  length,  behind  the  shock,  varies  as  Poo/Pgf 
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13 

For  constant  p|i|  it  has  been  found  that  the  interpolation  formula 

tj  =  3.  20/VI  +  0.  543W  is  accurate.  This  together  with  the  values  of  1^  in 

Table  A- 1  permits  Eq.  (55)  to  be  evaluated  for  constant  pp.  The  effect 

of  variable  pp  on  (rj^  -  1^)^  for  air,  has  been  determined  in  Reference  13. 

It  was  found  that  the  variable  pp  values  of  (n  -  I  )  for  air,  can  be  obtained 

u  00  0  0. 48 

from  the  constant  pp  values  by  multiplying  the  latter  by  (C  )  .It  may 
be  assumed  that  the  latter  correction  for  variable  pp  is  also  reasonable  for 
other  gases. 

Values  of  (26^/d)|m  have  been  computed  for  air  and  for  argon,  including  the 

effect  of  variable  pp.  It  was  assumed  that  (3  ?  Pj  in  Eq.  (55).  The  results 

are  given  in  Fig.  16.  It  is  seen  that,  for  both  argon  and  air, 

(26  /d )g  %  1,0.7,  and  0.4  for  M  1.2,  1.6,  and  2.  4,  respectively,  and 

u  m  8 

decreases  with  increase  of  M  . 

s 

The  excess  mass  flow  in  the  boundary  layer  was  previously  found  by  integrations 
in  which  it  was  assumed  that  the  boundary  layer  thickness  small  compared  with 
d  (e.  g.  ,  Eq.  6).  Since  most  of  this  mass  flow  occurs  in  the  portion  of  the 
boundary  layer  near  the  wall,  6^  is  a  conservative  criterion  for  the  size  of  the 
boundary  layer  relative  to  the  tube  radius.  For  laminar  boundary  layers,  the 
thin  boundary  layer  assumption  is  probably  valid  for  as  low  as  about  1.6. 
This  covers  the  range  of  Mg  of  practical  interest  in  low  pressure  shock  tubes. 
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VII.  SUMMARY  AND  CONCLUDING  REMARKS 


Test  time  and  flow  nonuniformity  in  low  pressure  shock  tubes  have  been 
investigated.  It  was  assumed  that  the  boundary  layer  was  laminar  and  was 
thin  relative  to  the  tube  diameter. 

In  the  first  portion  of  the  study,  the  asymptotic  flow  after  long  times  was 
considered.  Here,  the  shock  and  contact  surface  moved  with  constant  and 
equal  velocity.  The  flow  between  shock  and  contact  surface  was  steady,  in  a 
shock  fixed  co-ordinate  system,  and  the  separation  distance  was  found.  It 
was  necessary  to  treat  simultaneously  the  boundary  layer  development  and 
the  change  in  free  stream  conditions  external  to  the  boundary  layer.  A  local 
similarity  boundary  layer  solution  was  used  which  utilized  the  uniform  free 
stream  solutions  given  in  References  11  to  13.  The  local  similarity  solution 
should  be  accurate  except  for  Mg  near  1  (where  neglect  of  the  pressure 
gradient  on  the  boundary  layer  profile  may  become  important).  The  accuracy 
of  the  local  similarity  solution  can  be  established  by  more  accurate  solutions 
of  the  wall  boundary  layer  or,  perhaps,  by  applying  the  same  procedure  to 
boundary  layer  problems  where  the  solutions  are  known.  This  has  not  been 
attempted. 

The  asymptotic  separation  distance  yielded  values  of  p  which  were  consid¬ 
erably  larger  than  the  previous  analytical  estimates  of  Roshko  and  Hooker, 

particularly  for  moderate  M  .  (As  an  example,  consider  y  “  7/5,  <r  =  l,  and 

s 

pp  constant  across  the  boundary  layer.  Roshko  obtained  0.91  £  pj^  <  1.  52, 

for  1 . 29  <  M  <oo,  whereas  we  have  obtained  2. 22  <  p .  <_  1 . 81  for  the  same 

8  -2  A 
range  of  M  .  )  Test  time  is  proportional  to  p  .  Hence,  the  previous 

analytical  estimates  for  test  time  are  much  too  large,  particularly  for 

moderate  M^. 

The  experimental  test  time  data  of  Duff,  Roshko,  and  Hooker  indicated  that 
the  present  estimates  for  Pj,  are  still  somewhat  low.  For  argon,  they  are 
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too  low  by  about  35  percent  for  M  =  1 . 6;  20  percent  for  M  =  4,  and 

8  8 

10  percent  for  5  <  M  <9-  They  are  also  about  1 0  percent  low  for  5  <  M  <9 
8  8 

in  air.  The  discrepancy  between  theory  and  experiment  decreases  with 
increasing  Mg.  This  trend  is  to  be  expected  since  the  local  similarity  solu¬ 
tion  becomes  more  accurate  as  increases.  It  is  also  to  be  expected  that 
the  analytical  estimates  for  pj  should  be  lower  than  the  experimentally 
observed  values  since  mixing  and  diffusion  at  the  contact  surface  will  tend  to 
reduce  test  time  (i.  e.  ,  increase  p). 

In  the  second  portion  of  the  present  paper, consideration  was  not  restricted  to 

the  asymptotic  flow,  and  the  separation  distance  was  found  as  a  function  of 

time.  In  order  to  make  the  problem  tractable,  it  was  assumed  that  the  shock 

moved  with  uniform  velocity  and  that  the  flow  between  the  shock  and  the 

contact  surface  was  steady  (in  shock  fixed  co-ordinates).  Both  assumptions 

are  somewhat  in  error  but  should  give  at  least  qualitatively  correct  results. 

The  variation  of  separation  distance  with  time,  in  nondimensional  form,  was 

essentially  the  same  as  that  obtained  by  Roshko,  except  for  near  1.  Hence, 

the  main  difference  between  Roshko1  s  results  and  the  present  results  is  in  the 
2 

numerical  value  of  p  (which  is  used  to  nondimensionalize  both  separation 
distance  and  time). 

The  shock  Mach  number  varies  from  an  initial  value  of  M  .to  an  asymptotic 

s ,  1 

value  M  A  during  the  course  of  the  flow  in  a  low  pressure  shock  tube.  For 

strong  shocks,  Mg  ^  nearly  equals  Mg  ^  and  the  assumption  of  a  uniform 

shock  velocity  is  valid.  For  weaker  shocks,  the  question  might  arise  as  to 

what  value  of  Mg  should  be  used  to  correlate  experimental  test  time  data  or  to 

theoretically  predict  test  time.  When  correlating  experimental  data,  the 

locally  observed  value  of  Mg  should  be  used.  When  predicting  test  time,  the 

local  value  of  X  =  x  /Wf  should  first  be  evaluated  using  either  M  .  or 
s  m  °  s,  i 

M  A  (Eq.  50).  If  X  >  0(10),  the  local  flow  will  be  in  the  asymptotic  flow 
s ,  * * 

condition,  and  M  =  M  A.  If  X  <0(0.  1),  the  local  flow  is  only  slightly  per- 

S  S,  A 

turbed  from  the  ideal  flow,  and  M  =  M  For  X  =  0(1),  the  local  flow  is  in 

s  s,  1 

an  intermediate  condition  and  a  mean  value  between  M  .  and  M  A  would  be 

8  |  1  8 1  A 

appropriate. 
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APPENDIX 


BOUNDARY  LAYER  BEHIND  MOVING  SHOCK 


I 


The  results  given  in  References  11  and  12  are  summarized  with  regard  to  the 
displacement  thickness  of  the  boundary  layer  behind  a  shock  moving  with 
uniform  velocity.  It  is  assumed  that  y  and  <r  are  constant,  that  pp  is  constant 
across  the  boundary  layer,  and  that  the  free  stream  is  uniform  behind  the 
shock. 

These  results  are  then  put  in  a  form  applicable  to  the  local  similarity  solution 
discussed  in  the  body  of  the  report. 

UNIFORM  FLOW  BEHIND  SHOCK 

The  boundary  layer  is  indicated  in  Fig.  3.  The  free  stream  is  uniform  and  is 
the  same  as  that  directly  behind  the  shock.  Subscript  o  is  used  for  this  flow. 

In  References  11  and  12  a  similarity  variable  r|  is  employed,  defined  by 


dti 


and  a  stream  function  fQ(r|)  such  that 


f'  = 
o 


u/u 


e,  o 


£o(0)  =  o 


(A- 1 ) 


(A-2) 


The  boundary  layer  displacement  thickness  is  given  by 

which  is  negative  for  the  boundary  layer  behind  a  moving  shock. 
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In  terms  of  transformed  quantities, 
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+  I  ] 
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=  G 
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Equation  (A-4)  defines  (f  -  n)  ,  I  ,  and  G  .  These  quantities  are  of 

ao,  o  oo,  o  o 

fundamental  importance  in  determining  shock  tube  test  time  and  can  be 
evaluated  from  the  numerical  results  given  in  References  11  through  13. 
Typical  values  are  given  in  Table  A- 1  for  fluids  with  constant  specific  heats, 
constant  cr,  and  constant  pp.  Real  gas  results  for  air  are  given  in  Table  2. 

The  boundary  layer  parameters  in  Table  A-l  were  obtained  as  follows.  The 

quantity  (f  -  r\)  is  a  function  only  of  W  s  u  /u  and  is  tabulated  in 
oo i  o  w  e,  o  |  j 

References  11  and  12.  The  quantity  I  is  found  from 

n  7  oo,  o 
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where  M  =  u  /a  and 
e,  o  e,  o  e,  o 
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The  quantities  I  ,  I  ,  and  r  (0)  are  functions  of  cr,  as  well  as  W,  and  are 
^  r,  o  s,  o  o 

defined  and  evaluated  in  References  11  and  12  for  (r  =  0,  72.  However,  for 
<r  =  1,  these  quantities  become  (using  Ref.  11) 
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where  f£(0)  is  also  a  tabulated  function  of  W.  Normal  shock  relations  give 
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=  ZW"T"  l 


(A-7a) 
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e,  o 


2W  -  1 


(A- 7b) 


where  Z  =  (y  +  l)/(y  -  1)  In  Eq.  (A-7b),it  is  assumed  that  T  =  T  The 

2  V/  00 

shock  Mach  number  Mg  is  related  to  W  and  Z  by  Mg  =  W(Z  -  1)/(Z  -  W). 
Equations  (A-7a)  and  (A-7b),  combined  with  the  tabulated  results  in  Refer¬ 
ences  11  and  12,  yield  the  values  of  I  and  G  in  Table  A-l. 

7  oo,  o  o 

The  con&tant  p jjl  solutions  in  References  11  and  12  have  been  correlated  in 
Reference  12  by  interpolation  formulas.  The  equations 

(f  -  q)  =  1.  1 35(W  -  1)A/1  +  1. 022W  (A- 8a) 

00,  o  T 

f»(0)  =  -0.  489 (W  -  1)  yi  +  1. 665W  (A-8b) 


correlate  the  numerical  results  in  References  11  and  12  to  within  1  percent. 

It  is  also  of  interest  to  find  an  approximate  analytic  expression  for 
I  Q/(f  -  r))^  This  can  be  readily  done  for  a  =  1  by  using  Eqs.  (A- 5) 
through  (A-8).  We  find 


oo,  o 


=  [1  +  0.  431^(1  +  1.665WK1  +  1.022W]j/(ZW  -  1)  (A-9a) 
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For  strong  shocks*  Z  can  be  replaced  by  W  in  Eq.  (A-9b). 

BOUNDARY  LAYER  THEORY  FOR  LOCAL  SIMILARITY  SOLUTION 

The  boundary  layer  theory  for  use  in  the  local  similarity  theory  presented 

in  the  body  of  the  report  will  now  be  outlined.  For  simplicity,  it  is  assumed 

that  o*  a  1 .  The  free  stream  is  uniform,  with  a  velocity  uq»  and  the  wall  has 

a  velocity  u  (Fig.  4).  The  variable  V  =  u  /u  is  introduced.  As  before, 

7  w  e  e  e,  o 

W  =  u  /u  . 
w  e,  o 

Expressions  for  (f  -  r\)  and  fM(0)  are  found  by  replacing  W  by  W/V  in 
Eqs.  (A-8).  The  result  is 


(/  1. 13s(£  - 


(A-10a) 


f"(0)  =  -0.489^  -  l)yi  + 1.665  ^  (A- 10b) 

which  is  valid  for  all  <r.  For  <r  =  1,  Eq.  (A- 5)  becomes  (since  r(0)  =  1) 
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(A- 11) 


where 


1  = 
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I  = 


•  ($-0 


Bat  Me  =  „e/s  ■  Me  oV^Te  o/Te..d  Tw/Te  •  (T./T,  .XT.  0/T.).  Al.o, 

the  flow  in  the  free  stream  is  isentropic  so  that  p  /p^  _  s  (pmtpm  an^ 

.  e  e,  o  e  c,  u 

T  /T  =  (p  /p  )v.  From  Eqs.  (A-7),  (A- 10),  and  (A- 11),  it  can  then 
e  e  i  o  e  e ,  o 
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be  shown  that 
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It  also  follows  that 
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(A-l  3) 

To  integrate  Eq.  (16),  it  is  necessary  to  express  as  a  function  of  This  is 
done  as  follows.  For  isentropic  flow 
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which fwith  Eq.  (8),  gives 
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which  can  be  solved  by  iteration  to  find  M  as  a  function  of  7.  Also*  for 

e 

isentropic  flow 
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Finally 


V  = 


P  u 
e  e 

Tp  "u  r 

e  e  o 


e ,  o  _ 
P 


=  o 


(A- 15) 


(A-16) 


Equations  (A-13)  through  (A-16)  permit  to  be  found  as  a  function  of  T  and 
permit  the  integration  of  Eq.  (16). 
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